Measurement of the branching ratio of n° — > e + e using K L — > 37r° decays in flight 
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The branching ratio of the rare decay it — > e + e _ has been measured in E799-II, a rare kaon decay ex- 
periment using the KTeV detector at Fermilab. We observed 275 candidate it — > e + e~~ events, with an ex- 
pected background of 21.4 ± 6.2 events which includes the contribution from Dalitz decays. We measured 
BR(tt° -> e+e", (m e+e -/m^o) 2 > 0.95) = (6.09 ± 0.40 ± 0.24) x 10~ 8 , where the first error is statistical 
and the second systematic. This result is the first significant observation of the excess rate for this decay above 
the unitarity lower bound. 
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The decay 7r° — > e + e has received much experimental 
and theoretical attention since its branching ratio was first cal- 
culated by Drell in 1959 [0. Within the Standard Model, 
this decay proceeds predominantly through a two-photon in- 
termediate state, at a rate less than 1CP 7 that of ir° — ► 77. 
Relative to ir° — ► 77, ir° — ► e + e~ is suppressed by a helic- 
ity factor (2m e /m 7r o) 2 as well as by two orders of «em- 
The contribution to the decay from on-shell internal photons 
has been calculated exactly in QED [^], and forms the "uni- 
tarity bound," a lower limit on the branching ratio which is 
BR(7T° — > e + e~) > 4.69 x 10~ 8 , ignoring final-state radia- 
tive effects. The contribution from off-shell photons depends 
on the 7r°7*7* form factor, and is model-dependent. Recent 
vector meson dominance [||] and chiral perturbation theory 
calculations predict this contribution to be somewhat 
smaller than that from on-shell photons, giving a total branch- 
ing ratio of (6 — 9) x 1CP 8 , ignoring radiative corrections. 

Earlier experiments have produced conflicting measure- 
ments of the branching ratio for this mode. The earliest 
measurements of ir° — ► e + e~ were performed by a Geneva- 
Saclay group in 1978 [|| using 7r°'s produced by the decay 
K+ -> tt+tt in flight, and by a LAMPF group in 1983 
|0] using the charge-exchange process ir~p — > ir°n from a 
300 MeV/c pion beam. Both experiments favored a branch- 
ing ratio of ~ 2 x 10~ 7 , which would be hard to accomodate 
within the Standard Model. A 1989 search by the SINDRUM 
collaboration [H], using stopped ir~p — > 7r°rt, produced a 90% 
confidence level upper limit of 1.3 x 10~ 7 , excluding the cen- 
tral values of both previous measurements. In 1993, BNL 
E851 [§] and FNAL E799-I [Q observed the decay at the 
(5 — 10) x 10~ 8 level, near the Standard Model expectation. 
The BNL measurement used K + — > 7r + 7r° decays, while the 
FNAL experiment used Kl — > 37r° decays. 

In this Letter we present a new, precision measurement of 
BR(7T° — > e + e~) from E799-II, a rare Kl decay experiment 
which took data in 1997 at Fermilab. The 7r°'s were produced 
using Kl — > 3ir° decays in flight, where the other two 7r°'s 
in the event decayed to 77. The ir° — ► e + e~ events were 
normalized to Dalitz decays (ir° — ► e + e~j) with m e + e - > 
65 MeV/c 2 , which were collected and analyzed simultane- 
ously. High-m e + e - events were used in order to keep the 
charged track kinematic variables as similar as possible for 
the signal and normalization modes, and thus cancel many 
detector-related systematic errors. 

This technique, which was adapted from E799-I, has signif- 
icant advantages over those used in other measurements. The 
7T° 's were produced and decayed in vacuum, eliminating back- 
grounds and resolution smearing from decay products scatter- 
ing or converting in charge-exchange targets. The continuum 
process Kl — * 7r°7r°e + e~ has never been observed and does 
not pose a significant background. By contrast, the analogous 
processes in K + experiments (K + — > ir + e + e~) and charge- 
exchange experiments (n~p — > e + e~X) both produce large 
backgrounds to ir° — > e + e~. Reconstruction of the full kaon 
decay provides redundant kinematic constraints, eliminating 
all non-i^L — > 37r° backgrounds. 



The elements of the E799-II spectrometer (Fig. [T]) relevant 
to this measurement are described below. Two nearly-parallel 
neutral kaon beams were produced by 800 GeV protons strik- 
ing a 30 cm BeO target at a targeting angle of 4.8 mrad. Two 
neutral beams, each one up to 0.35 /isr, were defined by col- 
limators. A Pb absorber converted photons in the beam, and 
charged particles were removed by a series of sweeping mag- 
nets. A 65 m evacuated decay region ended at a Mylar-Kevlar 
vacuum window 159 m from the target. The beams in the de- 
cay region were composed mostly of neutrons and Kl, with 
small numbers of K$, A , A , 3°, 3°. These short-lived par- 
ticles tended to decay upstream. The Kl momentum ranged 
from - 20 to - 200 GeV/c. 

Charged particles were tracked using four drift chambers 
with two orthogonal views; a dipole magnet downstream of 
the second chamber provided a transverse momentum kick of 
205 MeV/c. Helium bags were placed between the chambers 
to reduce multiple scattering and photon conversions. 

Photon energy measurement and electron identification 
were performed using a 3100-block pure Csl electromag- 
netic calorimeter. The photon energy resolution was ~ 
1%, averaged over the energy range typical of 7r° — > e + e~ 
events (2—60 GeV). Immediately upstream of the Csl, two 
overlapping banks of scintillation counters (the "trigger ho- 
doscopes") provided fast signals for triggering on charged par- 
ticles. Downstream of the calorimeter, a 10 cm lead wall fol- 
lowed by a scintillator plane formed a hadron veto which re- 
jected at trigger level events with charged pions in the final 
state. An eleven-plane photon veto system, consisting of lead- 
scintillator counters throughout the decay region and spec- 
trometer, detected particles which left the fiducial volume. 

The trigger required at least 24 GeV of in-time energy in 
the Csl, and hits in the drift chambers and trigger hodoscopes 
consistent with at least two tracks. Events were rejected when 
more than 0.5 GeV was deposited in any photon veto counter, 
or more than the equivalent of 2.5 minimum ionizing parti- 
cles were detected in the hadron veto. A hardware processor 
JTT| ] required at least four energy clusters in the Csl, where a 
"cluster" was a set of contiguous blocks with at least 1 GeV 
deposited in each block. 

Offline, events with exactly two reconstructed tracks were 
selected and the tracks were required to form a common vertex 
inside the decay region. They also had to be electron candi- 
dates, defined as tracks which pointed to a Csl cluster whose 
energy E was within ±8% of the track momentum p. The re- 
constructed kaon energy was required to be at least 40 GeV, 
and each cluster in the Csl at least 1.5 GeV. These cuts re- 
duced the dependence of the result on Csl trigger thresholds. 
Clusters without tracks pointing to them were assumed to be 
photons. 7T° — > e + e~ candidates were required to have four 
photons; 7r° — > e + e~j candidates were required to have five. 

Photons were reconstructed assuming that each pair of pho- 
tons came from a ir° — > 77 decay. We calculated the distance 
Z12 = (ri2 / m„o )y/E~\E~2 of the decay from the Csl, where 
Ei is the energy of photon i and is the transverse separa- 
tion of photons i and j at the Csl. For 7r° — > e + e~ candidates, 
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there were three possible pairings of the photons into two 7r°'s. 
The pairing was chosen which minimized the \ 2 f° r the hy- 
pothesis that the two ir° decays occurred at the same position. 
For tt° — > e + e~7 candidates, there were fifteen pairings; the 
best was selected and the unpaired photon was assumed to 
have come from the 7r° 



e~7 decay. 

The photon four-momenta were calculated assuming the 
photons originated at the weighted average of the Z posi- 
tions of the two 7r°'s and the transverse position of the recon- 
structed two-electron vertex. This reconstruction method al- 
lowed kinematic quantities to be calculated in nearly the same 
way for the signal and normalization modes, thereby cancel- 
ing certain systematic errors. 

The basic reconstruction cuts described below were ap- 
plied, and the samples obtained were used to study acceptance 
and backgrounds. The total invariant mass ra^a^o e + e - was 
required to be within 50 MeV/c 2 of the Kl mass. The to- 
tal momentum transverse to the kaon direction was required 
to be less than 30 MeV/c. For the normalization sample, 



the Dalitz decay mass m e 



was required to be within 



30 MeV/c 2 of the 7r° mass, and the reconstructed electron pair 
mass m e + e - > 70 MeV/c 2 in order to avoid systematic errors 
from mass resolution smearing near the 65 MeV/c 2 cutoff. 

A detailed Monte Carlo (MC) simulation was used to es- 
timate acceptance for the signal and normalization modes, as 
well as the level of backgrounds in the samples. Both the 
signal and normalization MC were implemented with radia- 
tive corrections. The tt° — ► e + e~ MC used the 0(o!em) radi- 
ation model of Bergstrom [ |T2| | , and the 7r° — > e + e~7 MC 
used an 0(a| M ) calculation based on the work of Mikaelian 
and Smith [Q. 

At this stage, the sample in the ir° — > e + e~ signal re- 
gion (0.132 < m e + e - < 0.138 GeV/c 2 ) was background- 
dominated. Fig. |] shows the distribution of m e + e - for data 
and for the MC background predictions. The backgrounds, 
which all came from Kl — > 5tt° decays, were as follows. 
Very high-?7j e + e - Dalitz decays (tt° — > e + e~j) could be mis- 
reconstructed as 7T° — > e + e~ if the photon was not detected 
and m e + e - was reconstructed slightly high (by 1-10 MeV). 
Another type of background came from decays with four 
electrons in the final state. When one electron of each sign 
was soft, the spectrometer magnet could sweep them out of 
the fiducial volume, leaving only two reconstructible tracks. 
The four electrons could come from Kl — > 3ir° with multi- 
ple 7T° — > e + e~7 decays, from a 7r° — > e + e~e + e~ decay, or 
from photon conversions in the (3.55±0.17) x 10 -3 radiation 
length vacuum window assembly [ 14 1. 

The four-electron backgrounds fell into two categories. 
1) "Correctly paired" four-track backgrounds, where all 
four electrons came from the same tt°: these included 
i) 7T° — > e + e~e + e~ decays, ii) 7r° — > e + e~7 where the pho- 
ton from the Dalitz decay converted, iii) 3ir° — > 67 where 
two photons from the same ir° converted. In correctly paired 
four-track events, the reconstructed m e + e - was generally be- 
low m T o, and m T o^o e + e - was slightly below mic L . 2) "Mi- 
spaired" four-track backgrounds, where the four electrons 



came from different 7r°'s: These included i) events where two 
7r°'s decayed to e + e~j, ii) 7r° — > e + e~j events where a pho- 
ton from a different ir° converted, iii) 3ir° — > 67 events where 
two photons from different ir° decays converted. In these 
cases, because the two observed electrons did not come from 
the same 7r°, the m e + e - distribution was nearly flat. Because 
the four photons were not from two tt° — ► 77 decays, the Z 
position and the photon four-momenta were misreconstructed, 
giving a flat m 7T o^o e + e - distribution as well. 

Requiring the pairing x 2 to be below 4.5 removed 88% of 
the mispaired four-track background at the cost of 10% of the 
signal. Tightening the total mass cut to |m 7r o 7r o e + e -—m^ i I < 
10 MeV/c 2 removed a further 80% of the mispaired four-track 
background with negligible signal loss. 

The correctly paired four-track background could not be 
reduced significantly with pairing or kinematic cuts. About 
99% of these events, as well as 98% of the remaining mis- 
paired four-track background and 8% of signal, were removed 
by cutting events with evidence of extra in-time tracks in the 
second drift chamber. The last three cuts were applied to both 
the signal and normalization samples. 

After all cuts, the total background was dominated by high- 
mass Dalitz decays (18.1 ± 4.7 events). Smaller backgrounds 
came from correctly paired (2.8 ±1.1 events) and mispaired 
(0.5 ± 0.5 events) four-track modes. The errors on the four- 
track backgrounds are from MC statistics; the error on the 
Dalitz background reflects MC statistics and a 20% system- 
atic uncertainty in the MC prediction of the misreconstructed 
m e + e - tail. An (18 ± 5)% discrepancy was seen between 
the data and the MC prediction in the level of the low-m e + e - 
Dalitz background between 0.110 < m e + e - < 0.125 GeV/c 2 
(Fig. [3]). Perfect agreement in this region was not expected, as 
these events typically had an extra low-energy photon near the 
cluster energy threshold. Although the Dalitz decays which 
entered the signal sample had a much lower-energy photon 
and were therefore less sensitive to the modeling of the thresh- 
old in the MC, we have treated this discrepancy conservatively 
as an additional systematic error on the background. The final 
background estimate was therefore 21.4 ± 6.2 events. 

Radiative corrections to ir° — > e + e~ had a significant ef- 
fect on the acceptance. Internal bremsstrahlung can produce a 
e + e~7 final state with m e + e - < m^o, indistinguishable from 
the tree-level Dalitz decay. Following the convention of Ref. 
JlCQ, we imposed a cutoff (m e + e - /m w o) 2 > 0.95. We thus 
quote the branching ratio for this range only, after subtract- 
ing the small contribution from the Dalitz diagram. In this 
kinematic region, interference between the two processes is 
negligible [@. 

The normalization data set contained 650 264 events, with 
negligible backgrounds. The acceptance for Dalitz decays 
with m e + e - > 65 MeV/c 2 was 1.03%, for kaons which de- 
cayed between 90 and 160 m from the target and had momen- 
tum between 20 and 200 GeV. The acceptance for the signal 
mode was 2.52%. 

In the data, 275 7r° — > e + e _ candidate events were ob- 
served (Fig. ||). Subtracting the estimated background yielded 
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the total sample size of 253.6 ± 16.6 events (the error is sta- 
tistical only). 

The largest source of systematic error was the 2.7% uncer- 
tainty in the n° — > e + e~7 branching ratio (1.198 ± 0.032) x 
10~ 2 p"5]]. In addition, the m e + e - cutoff in the normalization 
Dalitz decays introduced a dependence of the acceptance on 
the Dalitz decay form factor. The MC used the form factor 
slope of 0.033 ± 0.003 measured by the CELLO collabora- 
tion ||], which gives the result that the m e + e - > 65 MeV/c 
region contains 3.19% of all Dalitz decays. The CELLO mea- 
surement used the reaction e + e~ — > 7r°e + e~ in a region of 
spacelike momentum transfer, extrapolating the slope to the 
kinematic region of the Dalitz decay assuming vector me- 
son dominance. The most recent direct measurement from 
the Dalitz decay is consistent but less precise [|l7|]. Using the 
CELLO form factor, the observed m e + e - distribution (Fig. |]) 
was consistent with the MC. The statistical precision of our fit 
was 0.007, which we have taken to be the uncertainty in the 
form factor; this translates into a 0.5% systematic error on our 
measurement of ir° — ► e + e~. 

The Dalitz decay branching ratio and the background un- 
certainty dominated the systematic error. Smaller accep- 
tance uncertainties included a 1.0% uncertainty in the effi- 
ciency of the pairing x 2 cut and a 1.2% systematic error as- 
signed to the efficiency of the m e + e - 7 cut in the normaliza- 
tion sample. These errors were determined from resolution 
studies using fully-reconstructed Dalitz decays. Adding all 
the systematic errors in quadrature, we obtained a total sys- 
tematic error of 4.0%. Our result for the branching ratio is 
BR(tt -> e + e- ,(m e+e - jm^f > 0.95) = (6.09±0.40± 
0.24) x 10~ 8 , where the first error is statistical and the second 
systematic. 

For comparison with the unitarity bound and with theo- 
retical models which neglect final-state radiation, we can in- 
vert the radiative corrections and extrapolate our result to the 
"lowest-order" rate (what the branching ratio would be in the 

absence of final-state radiation). This yields Fe+ |7 = 

(7.04 ± 0.46(stat) ± 0.28(syst)) x 10~ 8 , which is over four 
standard deviations above the unitarity bound. This result, 
which is in agreement with recent Standard Model predic- 
tions, represents the first statistically signficant observation of 
an excess above unitarity. 

E799-II expects to accumulate 2 to 4 times more data in 
a 1999 run, which will allow a further refinement of this 
measurement. This result should provide constraints for pre- 
dictions of related decay modes such as rj — > and 
Kl —> M + M [Pl~|5yi8[l, and we hope that future experiments 
will be able to test these predictions. 
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FIG. 1. Plan view of the KTeV spectrometer as configured for 
E799-II. The horizontal scale is compressed. 
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FIG. 2. Reconstructed m e + e ~ after basic reconstruction cuts. 
Left: background MC predictions. Right: E799-II data overlaid on 
sum of background predictions. (Backgrounds are normalized to the 
observed number of fully-reconstructed 7r° — * e + e~7 decays.) 
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FIG. 3. Reconstructed m e + e - for data, background MC, and sig- 
nal MC after all cuts. The vertical dotted lines indicate the defined 
signal region (0.132 < m e+e - < 0.138 GeV/c 2 ). 
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FIG. 4. Reconstructed m e + e - 



for normalization Dalitz decays. 
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